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Non-positive, Markovian semigroups are sometimes used to describe the time evolution 
of subsystems immersed in an external environment. A widely adopted prescription to 
avoid the appearance of negative probabilities is to eliminate from the admissible initial 
conditions those density matrices that would not remain positive by the action of the 
semigroup dynamics. Using a continuous variable model, we show that this procedure 
leads to physical inconsistencies when two subsystems are considered and their initial 
state is entangled. 



1. Introduction 

The time evolution of a subsytem in weak interaction with an external environment 
can be modelled in terms of Markovian (i.e. memoryless) dynamical semigroups; 
these are families {74} of one-parameter (=time) maps acting on the states of the 
subsystem, typically described by density matrices, and obeying the forward in 
time composition law, ^t°ls = lt+s, t,s > 0. Such a description is very general and 
have been successfully used to describe various effects in o pen system dynamics, in 
quantum optics, quantum chemistry and atomic physicspJ - *^ 

Nevertheless, the derivation of such time evolutions from the microscopic inter- 
action of the subsystem with the environment is often based on ad hoc approxi- 
mations, justified a posteriori with various physical considerations. As a result, the 
obtained time evolutions might not be fully consistent: in particular, such dynamics 
do not in general preserve the positivity of the reduced density matrix. Excep- 
tions to this general result are obtained through the use of rigorous mathematical 
treatmentspJ - ^ 

These inconsistencies are either dismissed as irrelevant fo r a ll practical 
purposes^l or cured by adopting further ad hoc prescriptionsP^HISl In the lat- 
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ter case, the general attitude is to restrict the action of the non-positive dynamical 
maps to a subset of all possible initial reduced density matrices, those for which 
the time-evolution remains positive. This is equivalent to a suitable selection of the 
initial conditions for the starting state of the subsystem, a procedure sometimes 
referred to as "slippage of the initial conditions" . On physical grounds, this effect is 
viewed as the consequence of the short-time correlations in the environment, that 
have not been properly taken into account in the derivation of the Markovian limit 
of the evolution, usually based on weak-coupling limit techniques that get rid of 
transient effects. 

In the following, we shall reexamine this widely used prescription to cure pos- 
sible inconsistencies produced by non-positive, Markovian evolutions, and point 
out further potential problems of this approach. Some results on finite-dimensional 
systems have been reported in Ref. [17, 18]. In this note we shall instead examine 
the case of simple infinite-dimensional systems. More specifically, we shall study 
the behaviour of two independent bosonic oscillators, one evolving with a dissi- 
pative Markovian semigroup, while the other with a standard unitary dynamics. 
We shall limit our considerations to the phenomenologically relevant set of Gaus- 
sian statei^W^Il anc i to dissipative dynamics that preserve this set, the so-called 
quasi-free semigroups EZl— E9] 

Although very simple, this settings actually corresponds to a specific physical 
situation in quantum opticil^H^, that of two coherent states travelling in optical 
fibers, of which only one is subjected to noise, producing dissipative phenomena, 
while the other gives rise only to unitary, standard birefringence effects. The consid- 
erations that follow may therefore have direct phenomcnological and experimental 
relevance. 

We shall explicitly show that when the initial state of the two oscillators is 
entangled, redefining the initial conditions to make positive the single system time- 
evolution is not enough to cure all possible inconsistencies of the two-system dy- 
namics. Therefore, in order to have a physically acceptable time evolution for the 
two subsystems in presence of entanglement, the above mentioned procedure of 
restricting initial conditions should be tested against the presence of non-classical 
correlations. 

2. Single Mode Dissipative Dynamics 

We shall first study the dynamics of a single oscillator in weak interaction with an 
external environment, so t hat its reduced time evolution can be well represented 
by a Markovian semigroupPSl The states of this system will be represented by a 
density matrix p, i.e. by a positive Hermitian operator, with unit trace, acting on a 
bosonic Hilbert space TL. It can be identified with a Fock space, generated from the 
vacuum state through the action of polynomials in the creation and annihilation 
a operators; these operators obey the standard bosonic oscillator algebra: [a, a'] = 1, 
[aT, a'] = [a, a] = 0. 
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Among all densi ty matr ices p, the so-called quasi-free or Gaussian states are 
of particular interesl(12H221 as mentioned before, they can be easily produced in 
experiments in quantum optics. These states are defined by the property that the 
expectations (W(z)} = tr\W(z) p] of the Weyl operators W(z) = e za+z ° f are in 
Gaussian form, i. e. they are exponentials of a quadratic form in the complex variable 
z and its conjugate z*. 

Let us collect annihilation and creation operators in the column vector a and 
its hermitian conjugate row vector = (a* a). Then Gaussian states are fully 
characterized by the 2x2 matrix of bilinear correlations {p,,v = 1, 2)f23H2ll 

o„.fc..i)-*Mi,]-($> A)" (1) 

From its definition, it follows that the matrix G in ((T|) is non-negative: G > 0, so 
that (3 > 1 and (3(j3 — 1) > |a| 2 ; one can show that this is also a sufficient condition 
for the state p to represent a physical stateP^For simplicity, we henceforth assume 
(at) = (a) = 0; this condition can be easily released, by taking (W(z)) to include 
exponentials of linear terms in z and z*. 

As mentioned in the introductory remarks, our analysis is based on the assump- 
tion that the time evolution of a single oscillator immersed in a bath be Markovian 
and given by a semigroup, i.e. by a trace preserving, one parameter family of linear 
maps, acting on the set of density matrices p representing the oscillator states. These 
maps are generated by master equations of the following Kossakowski-Lindblad 
fonnJU-E 

Ml = -i[H,p(t)]+L[p(t)}, (2) 

with H an effective Hamiltonian; instead, L is a linear map representing the dis- 
sipative piece: it describes the effects of noise induced by the presence of the en- 
vironment. The semigroup evolutions that preserve the set of Gaussian states are 
generated by equation in form © with H and L that are quadratic in the creation 
and annihilation operators} 



H = a} , uj > 

L \p] = v{[ap, a f ] + [a,/9a f ]) + a([a^p, a] + [a\ pa]) 

-X*[a,[a,p}} -\[a*,[<J,p]] ; (3) 

the coefficients 77, a and A encode the physical properties of the environment and 
can be expressed in terms of the Fourier transform of the correlation functions in 
the bathffHH 

By using the two-dimensional vectors a and a^ introduced before, the dissipative 
term in ([3]) can be recast in the following compact form: 
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where the bath coefficients are now embedded in the 2x2 Kossakowski matrix: 

I, A* 



C =U „)■ < 5 > 

The dissipative parameters in ([5|) are not completely arbitrary. First of all, 77 and 
a need to be real in order to comply with the hermiticity preserving requirement 
of the generated semigroup. In addition, the request that the finite-time evolution 
7t generated by a master equation of the form [[2]). ([3| map physical states into 
physical states, i.e. preserve the positivity of the density matrix p(t) for all times, 
gives rise to further constraints. 

The simplest way to guarantee the positivity of the dynamics 74 is to require 
the Kossakowski matrix in ([5|) to be non-negative 

C > . (6) 

Indeed, this condition is equivalent to the requirement of complete positivity, a 
stronger property than simple positivity, that guarantees the physical consistency 
of the evolution 74 in all possible situations LiH^I Dynamics of such type, obtained 
from the equation @, with Hamiltonian as in §3§ and dissipative term as in 
(0, are known in the literature as quasi- free quantum dynamical semigroups l^ZMSSl 
as already observed, they are characterized by the property of transforming the set 
of Gaussian states into itself. 

Nevertheless, in phenomenological applications, the condition of complete pos- 

ToTl 

itivity is often dismissed as unnecessary^ 1 , since in many naive derivations of the 
equation ([2]) generating the reduced dynamics 74, the hierarchy among t he d issi p a- 
tive parameters imposed by the condition ([6]) appear not to be satisfied!12H16| i n 
such cases, it is easy to show that the resulting dynamics jt is not positive, hence 
physically unacceptable. 

In order to explicitly show this, it is convenient to consider the time evolu- 
tion of the correlations in (fTJ), as induced by the dynamics of the density matrix 
p(t), through the duality map: Trfa^ aj p(t)\ = Tr[(a M a£)(t) p\ . One explicitly 

&G(t) = At-G(t) + G(t)-A + B, (7) 

where the 2x2 matrices A and B contain the dependence on the hamiltonian oj 
and dissipative parameters 77, a and A: 



2 \ o — Tj — 2iuj J \ —A a 

The positivity condition G(t) > needs to be preserved for all times in order for 
the the evolution j t '■ p(Q) — > p[t) to be positive and thus physically acceptable. In 
order to prove that this is not the case when ([6|) is not satisfied, i.e. when r\a < |A| 2 , 
it is sufficient to prove that G(t') < 0, for a certain t' . Let us assume that indeed 
this is the case; since we start with a physical state, i.e. G(0) > 0, there will be 
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a time t" < t' for which G(t") has a null eigenvalue. Since we are dealing with 
evolutions of semigroup form, without loss of generality, we can set t" = 0. 

Let us now consider the average Q(t) = ( 1 J r |G(i)|\I r ), where the two-component 
vector is the eigenvector of G(0) corresponding to its null eigenvalue. Recalling 
([!]), one easily finds (^| = (—a*, f3), with the condition j3(f3—l) = |a| 2 , that assures 
the vanishing of the determinant of G(0). Clearly, since Q(0) = 0, G(t) will start 
acquiring a negative eigenvalue as soon as: 



Since r), u > and (3 > 1, in order to satisfy (fT0|) one needs to choose the phase of 
a in such a way 7Ze(a\) be large and negative, i.e. TZe(aX) = —\a\ |A|. Then, one 
can check that for r\ > cr, the condition (| 1 Ojl is always satisfied; this remains true 
also in the case r\ < er, provided |A| < a. As a consequence, when the Kossakowski 
matrix in ([5|) is negative, the evolution generated by ([2]), ((3J results non-positive, 
and therefore physically unacceptable. 

As mentioned in the Introduction, in order to cure this pathology and continue 
to use semigroups of non-po sitive m aps to model reduced dynamics, an ad hoc pre- 
scription has been proposed Q3HEI restrict the possible initial states p(0) to those 
for which p(t) = 7 t [p(0)], t > 0, as generated by j2|), j3|), remains a physical Gaus- 
sian state. The argument is that any Markov approximation neglects a certain initial 
span of time, namely a transient, during which no semigroup time-evolution is pos- 
sible due to unavoidable memory effects. The transient dynamics can be effectively 
depicted as a slippage operator that, out of all possible initial states, selects those 
which do not conflict with the Markov time-evolution when it sets in. As we shall 
see in the next section, the slippage-argument may cure the positivity-preserving 
problem, but appears to be inconclusive when dealing with bi- and multi-partite 
open quantum systems because of the existence of entangled states. 

3. Two-Mode Dynamics and Entanglement 

We shall now extend the treatment discussed so far to the case of two oscillators, one 
of which still immersed in an external environment and therefore evolving with the 
dissipative, semigroup dynamics generated by ([2]), |3]), while the other undergoes 
a unitary evolution. As mentioned in the Introduction, such a system can actually 
describe physically implementable situations in quantum optics, e.g. that of polar- 
ized, coherent states propagating along two fibers, one which is subjected to noise, 
while the other is not. Having a physically consistent theoretical treatment of the 
behavior of such setups is therefore of great phenomenological relevance. 

In this case of bipartite, two-mode settings, Gaussian states p are still charac- 
terized by having expectation (W(z\, Z2)) — tr[VF(zi, z-i) p\ of the Weyl operators 



&Q(0) < . 
Using (JT]) and JSJ, this condition reduces to: 



(9) 



r?(/3 - 1) +(7/3 + 21le{a\) < . 



(10) 
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W(zi,z 2 ) — e Ziai + Z2Cl2+z i a i +z 2 a 2 in Gaussian form. These states are again fully 
characterized by the correlation matrix G, which is now 4x4, and still defined as 
in ([1]). The column vector a and its conjugate a* are now four-dimensional, with 
components a^ = {a\,a\,a\,a2), p = 1,2,3,4, where a\, <2j, i = 1,2, are single 
mode creation and annihilation operators, obeying the standard oscillator algebra: 
[a,i, a]] = 6ij, [a], ajj] = = [a», a,-]. 

The time evolution is again a Markovian semigroup generated by a master equa- 
tion of the form @. The Hamiltonian is now the sum of two pieces: 

H = -ui{a\, a%} + -u; 2 {4, a 2 } , (11) 

while the dissipative contribution can be written as in ((4]) , with a 4 x 4 Kossakowski 
matrix with null entries, except in the upper left corner where it coincides with ([5]). 
Such compound dynamics T t = "ft ® Ut , being represented by the tensor product of 
the dissipative evolution "f t discussed in the previous section for the first oscillator 
mode and of a unitary one Ut for the second mode, preserves the Gaussian form of 
the states p. 

As in the single mode case, the time evolution equation can be equivalently 
formulated in terms of the correlation matrix; explicitly: 

d t G(t)=A^-G(t) + G(t)-A + B, (12) 

where A is diagonal with entries h{a — rj + 2uo\, a — rj — 2iu\, 2iuj2, —2iui2), while 
the inhomogeneous term is given by 

-('!!)■ 

with representing the 2x2 zero matrix and B is as in (JSj). The matrix G(t) can 
similarly be written in block form: 

a, (A wm\ 

(14) 



G(t) _(G 1 (t) W(t)\ 



where Gi(t) and G 2 (t) are 2x2 correlation matrices corresponding to the reduced 
single-mode Gaussian states obtained by tracing over the degrees of freedom of the 
second, Tr 2 [p(t)], respectively first, Tii[p(t)], oscillator. As such, they must satisfy 
the positivity condition: G»(i) > 0, i = 1,2, for all times. 

Therefore, to avoid appearance of negative probabilities when the single mode 
dynamics "ft is not positive, i.e. when r\a < |A| 2 , in line with the "slippage" prescrip- 
tion, one has to restrict the choice of the initial condition p(0) to those Gaussian 
density matrices for which the reduced states Tr2[/o(0)] gives rise to positive correla- 
tion matrices, Gi(i) > 0. Although this prescription clearly works for initial states 
in factorized form, it fails in general when p(0) is entangled. 

In order to make the incompatibility between slippage and entanglement appar- 
ent, let us consider an initial correlation matrix of the form: 

/ c< r~ ,1 / 2 r~ ,1 / 2 \ 
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with the conditions Gj > 0, i = 1, 2; note that the matrix in (fT5|) satisfy the physical 
condition G(0) > by construction. Further, the 2x2 correlation submatrices 
are chosen in such a way that their positivity is preserved under the action of 
the compound dynamics T t . Let G2 be invertible. Then, the 4x4 matrix in (| 1 5|) 
possesses a zero eigenvalue, whose associated eigenvector can be written as 



i*> = 



-g 2 - 1/2 g1 /2 |v>). 



(16) 



where \ip) is an arbitrary two-dimensional vector. Let us then consider the average 
Q(t) = I G Since Q(0) = 0, the matrix G(t) will start developing negative 
eigenvalues as soon as 



a f Q(o)<o 



(17) 



Explicitly, d t Q(0) = = (tp\B\ijj); therefore, the condition (fTTJ) can always 

be satisfied by choosing for the two- vector the negative eigenvector of B. 

The physical inconsistency of dynamics for which the Kossakowski matrix in (0 
is not positive is now apparent: the slippage prescription make them positive, but 
not completely positive; therefore, when the subsystem is bipartite, one can always 
find an entangled Gaussian state that is mapped out from the set of physical states 
by the compound dynamics T t . 

That the initial state p(0) giving rise to the correlations in (fTS)) is en tangled 
can be directly verified using the operation of partial transpositional^, i.e. of 
transposition involving only one of the two oscillators, say the second one. It can 
be most easily implemented on symmetric correlation matrices: 



V 



4} 



1. 



-Tr 



(18) 



where it amounts to the exchange a 2 <-> a 2 - This transformation clearly maps states 
into states for separable ones, but not in general for correlated ones: it thus provides 
a sufficient criterion for bipartite entanglement,which turns out to be also sufficient 
for two-mode states, as in in the present case!^ 

For the case at hand, one finds the the operation of partial transposition with 
respect to the second system results in the following transformation of the 4x4 
symmetric covariance: 



V — » V = T • V • T , 





0"! 



(19) 



where 1 is the 2x2 unit matrix, while ct\ is the first Pauli matrix. Equivalently, for 
G(0) in Q3J), one finds: 



G(0) -> G(0) 



Gi 

l/2„l/2 



,1/2q1/2 

G 2 



(20) 
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where G2 indicates transposition of G2. By construction, the initial density matrix 
p(0) corresponds to a physical state, i.e. G(0) > 0; if one further finds: 

G(0) < , (21) 

then the state p(0) is surely entangled. One can show that (|2"Tj) is equivalent to the 
following condition involving 2x2 submatrices of G(0)P§1 

Gi - G{ /2 Gl /2 [G 2 ] _1 Gl /2 G{ /2 < , (22) 
which in turn reduces to: 

[G^-ax [G 2 ] _1 CTl <0 . (23) 

Since G2 can always be written as in {l}, with /3(/3 — 1) > \a\ 2 , by direct com- 
putation, one immediately finds that the combination in the l.h.s. of (|23[) always 
possesses a negative eigenvalue, and therefore that p(0) is indeed an entangled state. 



4. Discussion 

The dynamics of a subsystem in weak interaction with an external environment can 
be described in terms of semigroups of linear maps jt generated by a Markovian 
master equation. Although rigorously proven only in special cases, this result is 
nevertheless believed to hold in general, since memory effects should disappear from 
the reduced subsytem dynamics as soon as all correlations in the environment have 
died out. 

This mathematical effective description of reduced dissipative dynamics needs to 
satisfy further physical constraints. In the first place, it must preserve the positivity 
of any initial density matrix: only in this case, the time-evolution turns out to be 
consistent with the interpretation of the state-eigenvalues as probabilities. In the 
second place, one has also to care of possible couplings with another system, not 
subjected to noise, and therefore to guarantee the positivity-preserving character 
also of the semigroup of maps of the form T t — 7* ® Ut, as studied in the previous 
section; this is only assured in all possible situations if 7t is completely positive. 

Nevertheless, in many applications, complete positivity is often dismissed as 
unnecessary. This attitude is based on the fact that this property guarantees full 
physical consistency of the reduced dynamics at the price of order relations among 
the parameters describing dissipation (c/. ((6|)), which may appear to have little 
physical motivations. Further, most phenomenological derivation of reduced dissi- 
pative dynamics lead to semigroup of linear transformation that are not completely 
positive. 

Unfortunately, such derivations often lead to time evolutions that are not even 
positive. To avoid inconsistencies, one usually restricts the possible initial states to 
those for which "f t remains positive (the so-called "slippage of initial conditions"). 
This prescription works also in the case of the evolution T t for two subsystems, 
provided the initial state is in separable form: p(0) = ^2iPiPi ® Pi \ Pi > 0; 
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J2iPi = 1> where p- 1 ^ and p^p are admissible states for the first and second subsys- 
tems, respectively. 

On the contrary, as shown in the previous section, when the initial state p(0) is 
not in factorizcd form, the evolved matrix p(t) — T t [p(0)] fails to be positive for all 
times. In keeping with the same attitude adopted for a single subsystem dynamics 
7t, to cure this additional inconsistency one has to further restrict the domain of 
applicability of T t . However, this is again a temporary solution; indeed, the whole 
discussion needs be repeated when three or more subsystems are considered. 

These considerations seems to suggest that there is an intrinsic incompatibility 
between the existence of entangled states and the slippage prescription adopted to 
cure the inconsistencies that non-completely positive, or even non-positive, reduced 
dynamics might produce. Note that this conclusion can not be dismissed as being 
purely academic; on the contrary, it seems to have a direct experimental relevance: 
as mentioned before, the scenario studied in the previous sections might have an 
actual realization in quantum optics. From this perspective, the widely used cure 
of redefining the initial conditions in case of non-positive Markovian dynamics does 
not appear to be completely satisfactory. 
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